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Abstract 

The pion structure function is investigated in a simple model, where pion and 
constituent quark fields are coupled through the simplest pseudoscalar coupling. 
The imaginary part of the forward 7*7r — > 7*7r scattering amplitude is evaluated 
and related to the structure functions. It is shown that the introduction of non- 
perturbative effects, linked to the size of the pion and preserving gauge invariance, 
allows a connection with the quark distribution. It is predicted that higher-twist 
terms become negligible for Q 2 larger than ~2 GeV 2 and that quarks in the pion 
have a momentum fraction smaller than in the proton case. 
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1 Introduction 

Deep-inelastic scattering (DIS) experiments provide us with a wealth of infor- 
mation about the structure of hadrons, usually cast in the form of structure 
functions. These data are only partly understood in the framework of quan- 
tum chromodynamics (QCD). Indeed, perturbative QCD is consistent with 
the Q 2 evolution of the structure functions at sufficiently high Q 2 and x [1]. 
However, it is neither able to predict the structure functions themselves, as 
the latter are supposed to result from non-perturbative effects, among which 
confinement and spontaneous chiral symmetry breaking, nor the magnitude 
of the initial value Qq from which the Q 2 dependence can be evaluated . The 
interest has recently widened to off-diagonal parton distributions [2-5], which 
potentially offer to reach complementary information, especially about parton 
correlations. 

Phenomenological quark models, which possess some non-perturbative as- 
pects and which are rather successful in reproducing low-energy properties 
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of hadrons, are expected to help us to understand the connection between 
DIS data and non-pertubative aspects. Pions and other low-mass mesons are 
the simplest systems for which effective models exist that incorporate, in some 
simplified way, such special QCD features such as spontaneous chiral symme- 
try breaking and anomalies. 

Our original plan was to investigate the properties of off-diagonal parton dis- 
tributions in a simple model for light mesons. However, we realized that the 
theoretical investigation of the diagonal distributions using phenomenological 
quark models is far from being settled. There have been several theoretical 
investigations along these lines in recent years [6-9]. They rely on the as- 
sumption that distributions evaluated in leading-twist approximation at small 
Q 2 , where these models apply, can serve as input in the DGLAP evolution 
equations [10] to generate parton distributions that are directly comparable 
with experimental data at large Q 2 [11-15]. This is a rather tricky point as it 
is not clear that a good approximation at low Q 2 can be evolved by perturba- 
tive equations to large Q 2 , as it is not sure that, in this regime, the forward 
7*7r — > 7*7r amplitude can be parametrised in terms of parton distributions 
(despite the existence of low-Q 2 parametrisations [14,15]). 

Some of the theoretical works mentioned above have been carried out in the 
framework of the Nambu and Jona-Lasinio (NJL) model [16]. This model is 
indeed quite successful in reproducing low-energy phenomenology and embod- 
ies chiral symmetry breaking, which is believed to be crucial for this success. 
However, investigations of the structure functions have given rather different 
results [7-9]. This situation originates from the fact that the NJL model needs 
to be regularized and that different regularizations yield different results. 

We show in this paper that these complications may in fact be avoided. We 
investigate the simplest model of a pion, in which the qqjr vertex is repre- 
sented by the simplest pseudoscalar coupling. Of course, the quark fields and 
the coupling constant should reflect in some way the properties of actual pi- 
ons. One possible guide is provided by the large N c limit of the NJL model. 
The latter is equivalent, in the chirally broken phase, to a cr-model with mas- 
sive constituent quarks and an effective pion-quark-quark coupling constant. 
Our starting point for the investigation of the parton distributions of the pion 
consists in a Lagrangian that includes massive pion and massive quark fields 
interacting through the simplest pseudoscalar vertex. We show below that 
there is no need to regularize this model as the imaginary part of the forward 
7*7r — > 7*7r scattering amplitude is finite. However, the interaction between 
the constituent quarks should vanish when their relative momentum is large 
enough. It is possible to cope with this requirement by imposing a finite mo- 
mentum cut-off, which mimics the effect of a pion wave function. Such a step 
leads to the appearance of a straightforward relation between the 7* — n cross 
sections and the quark distributions at high enough Q 2 , as we shall show 
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Fig. 1. Simplest diagrams contributing to the imaginary part of the forward ampli- 
tude for the scattering of a virtual photon by a neutral pion. Upper (lower) diagrams 
are referred to as box (crossed) diagrams. Dashed lines represent the discontinuity 
of the amplitudes or their imaginary parts. 



below. 



2 The model 



We consider an isospin triplet pion field n = (7T + ,7r°,7r ) interacting with 
quark fields tp through the Lagrangian density 

£int = ig{i> rj 5 ip). 7T, (1) 



where f is the isospin vector operator. It is our purpose to calculate the imag- 
inary part of the forward elastic 7* — it scattering amplitude, or equivalently 
the total cross-section, in the simplest approximation and extract from it the 
pion structure functions W\ and Wi- Finally, we want to see whether they are 
reducible to quark distributions. We give below explicit results for the neutral 
pion. 
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The relevant diagrams contributing to the imaginary part of the forward elastic 
scattering amplitude, up to first order in a, the fine structure constant, and 
to second order in g, are given in Fig. 1. We define the pion 4-momentum 
as p, the photon 4-momentum as q, the (u and d) quark constituent mass as 
m q and we use p 2 = m 2 , q 2 = —Q 2 , p.q — v, x — Q 2 /{2v). This leads to 
s = m 2 + Q 2 {l/x — 1). The imaginary part of the amplitudes can be written, 
using Cutkosky rules, as 

ImT^ = -2Cg 2 J d 4 l{t] ll/ D 1 D 4 D 1 D 2 }, (2) 
ImT^ = -2Cg 2 J d 4 l{t 2 tl/ D 1 D 4 D 3 D 2 }, (3) 



with the fermionic traces 



(7./ + m g )7 5 (7.(/ + p) + m g )7 5 (7./ + m q )\ , (4) 

*)L = Tr hi*(l-( l - Q) + ^9)75(7- (l + p-q)+ m q )lv 

(-f.(l+p)+ m 9 )7 5 (7.Z + m q )\, (5) 



and the fermion propagators (or their contribution to the cuts) 



J D 1 - 1 = f - m 2 q , D 2 = 2n5 [(I + p) 2 - m 2 q 

£,-1 = (J + p - q )2 _ m 2 ? £, 4 = 2?r5 ( (/ _ g)2 _ m 2) 



(6) 



The constant C = 5e 2 /(3847r 6 ) accounts for flavour, charge and loop momen- 
tum integration factors in the particular case (neutral pion) under consid- 
eration. Diagrams (b) and (d) have the same contributions as (a) and (c), 
respectively. The second ones are obtained from the first ones by the substi- 
tution I — > q — p — I in the integrals. Using T 1 = T a + T b , T 2 = T c + T d and 
T = T 1 + T 2 , one can rewrite 

ImT^ = -2Cg 2 J dH{t\ v D x D±D x D 2 + (l^ q - p -l)}, (7) 
ImT 2 , = -2Cg 2 J dH{t\ v D x D,D z D 2 + (I - q - p - I)}, (8) 



which makes gauge invariance (q^lmT^ = ImT^q" = 0) explicit by inspection 
of the integrand. It is worth emphasizing that the sum of the four diagrams is 
gauge invariant, but none of them alone has this property. 
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We perform the integrations in Eqs. 7 and 8 using Sudakov variables. We 
define / = — £p + r)q + l t , with l t perpendicular to both p and q. A little algebra 
shows that 



dH 5((l + pf - ml)5{{l - qf - m 2 q ) = f g , (9) 



with the change of variable 



(I 2 - m 2 - 2mV)(2x - 1) , x 



and with r\ and I 2 related to I 2 through the relations 



_ -I 2 + m 2 + ml (2x - 1) m 2 n (-l 2 + m 2 ) 
2i/ + Am 2 x v(2v + Am 2 x) ' 

(1 — Am 2 x(l 2 — mix) + (/ 2 — m 2 ) 2 ) 

2z/ + 4m 2 x 
m 2 (/ 2 — 4m 2 i(i — 1) — m 2 — mix) 

2v + 4m 2 a; 
m 2 (/ 2 — m 2 ) 2 
4z/ 2 + 8z/m 2 x 



The integration over the azimuthal angle is trivial and one is left only 
with the integration over the variable r = —I 2 . The integration bounds on 
r come from the positivity of the energies of the on-shell intermediate states 
(see Fig. 1) and from the space-like definition of l t (yielding l 2 < 0), which 
introduces the most stringent constraints. The latter take the following forms: 



t <2u + 



((1 — x)xml — (2 — x)m 2 q ) 
1 — x 

(m 2 - x 2 ml)(m 2 q - (1 - x) 2 m 2 ) 
2v{l-x) 2 ' 



(13) 



and 



t > 



(ml - (1 - x)ml) (ml - x 2 ml)(m 2 - (1 - xfml) 



+ 



x 



2v(l-xf 



(14) 



where only the first terms in the expansion of the r.h.s. in \ jv are given (the 
full expressions are used in numerical evaluations). 
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The structure functions W\ and W 2 being defined by 



= — ImT^ 
In 



they can be calculated through the contractions 

-6xWi + (u + 2m 2 7r x)W 2 



2x 



fJLV 



-2x(v + 2m 2 n x)W 1 + (is + 2m\x) 2 W 2 
Ax 2 



(15) 

(16) 
(17) 



We also checked that we get the same results for W\ and W 2 by computing 
directly the 7*7r — > qq transverse and longitudinal cross-sections. 

Although we give below results for the exact expressions, it is worth giving 
the large-Q 2 limit: 



W 1 



5g 2 

24vr 2 



In 



'2(1 - x y 

M 2 



ml 



(18) 



with M 2 = m 2 q 



m 2 x(l — x). 



At this point, the pion cannot be interpreted as a collection of partons with 
a probability distribution. Indeed, the crossed diagrams are not suppressed 
by a power of Q 2 compared to the box diagrams and therefore do not allow 
such an interpretation. This is manifest when their respective contributions 
are written down in the combined small-m^, large-Q 2 limit: 



box 



24vr 2 



24tt 2 ' 



(19) 



A reason for the appearance of this undesirable feature is that we have not yet 
imposed the fact that the pion has a finite size. The simplest way to do this is 
to require that the square of the relative four-momentum of the quarks inside 
the pion is limited to a maximum value A 2 . The former quantity is given by 

Ox = ip + 2lf = -2t + 2m 2 , - m 2 , (20) 

for vertices like those in diagram Fig. l.(a), and by 

2 = (p + 2(1 + p - q)) 2 = 2t + 6m 2 - m 2 - 4z/, (21) 
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Fig. 2. Values of the coupling constant g(Q 2 ,A) which fulfill sum rule (Eq. 25), for 
the values of the parameters indicated at the top. 



for vertices like those in diagram Fig. l.(b). We either require |Oi| < A 2 or 
| O2 J < A 2 for all diagrams, which excludes the interval [A 2 , 2is — A 2 /2] from 
the r integral. 

I Oi| and 1 02 1 cannot be small simultaneously. The crossed diagrams have their 
main contribution for 0\ ~ O2, and are thus suppressed by a power A 2 /Q 2 
when the cut is imposed. The box diagrams have a leading contribution for 
1 0i I or 1 2 1 small, and are not power suppressed by the cut. 

Physically, this happens because, for the crossed diagrams, the momentum 
transfer suffered by the quark (antiquark) has to be re-emitted by the an- 
tiquark (quark), which is impossible when this momentum transfer becomes 
too large. For the box diagrams, the transferred momentum is taken and re- 
leased by the same quark (or antiquark) and there is no suppression. Note 
that t = —t (for the 7*7r — > qq process) is a physical observable quantity. This 
guarantees that our implementation of the cut-off is gauge-invariant. With our 
cut-off, the crossed diagrams appear now as higher twists: they are suppressed 
at least as 1/Q 2 and the box diagram contribution can now be interpreted in 
terms of parton distributions. Indeed, we checked numerically that, for typical 
values of A (see Fig. 2), the ratio between the respective magnitudes of the 
crossed and box diagrams is roughly ~0.06/<5 2 GeV~ 2 , for x—0.1; without 
cut-off, this ratio is a sizable fraction of unity for Q 2 less than ~ 1000 GeV 2 . 
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Fig. 3. Asymptotic value (for large Q 2 ) of the coupling constant fulfilling sum rule 
(Eq. 25), as a function of the cut-off parameter A. 

3 Results 



The structure functions can now be related to the (valence) quark distribu- 
tions: 



4 1 

F 1 = W 1 = —(u v (x) + u v (x)) + —(d v (x) + d v (x)), 

lo lo 



F 2 = uW 2 = 2xF 1 . 



(22) 



(23) 



We stress that the last relation, known as the Callan-Gross relation, comes 
out of our calculation. This does indicate that our approximations have been 
done consistently^ 

So far, we have only described the model. In order to make predictions, we 
need to fix its parameters, namely A, m q and g. The latter can be thought 
of as the normalisation of the quark wave function, and we determine it by 
imposing that there are only two constituent quarks in the pion. In our model, 
the valence quark distributions are equal 



u v (x) = u v (x) = d v (x) = d v (x) = v(x). 



(24) 



For charged pions, the additional diagrams implying a direct coupling of the 
virtual photon to the pion are suppressed by a factor 1/s and the leading-twist 
results are the same, except for the charge coefficients entering Eq. 22. 
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Fig. 4. Momentum fraction of the quarks inside the neutral pion (Eq. 26), as a 
function of Q 2 , for the values of the parameters indicated at the top of the figure. 

The condition v(x)dx = 1/2 then gives us 



J F^dx 



18' 



(25) 



As Fx is a function, not only of m q and A, but also of Q 2 , this gives us a 
coupling constant that evolves with Q 2 . The resulting values of g are shown 
in Fig. 2: for a finite cut-off A, the cross-section at fixed g would grow with 
energy, until the pion reaches its maximum allowed size, in which case the 
cross section would remain constant. If we impose relation (25), this means 
that g{Q 2 ) will first decrease until the cut-off makes it reach a plateau value 
for Q 2 3> A 2 (in practice, the plateau value is reached around Q 2 2A 2 ). The 
plateau value depends on m q /A (and m n /A) and is shown in Fig. 3. 

To constrain further our parameters, we can try to use the momentum sum 
rule 



f 18 f 

2 (x) = 4 J xv(x)dx = — J F 2 (x)dx. 



(26) 



In the parton model, this integral should be equal to 1 as Q 2 — > oo, as we 
do not have gluons in the model. Indeed, there are two limits of our model 
that fulfill this condition automatically. First of all, if we do not impose that 
the pion has a finite size, then asymptotically we get 2 (x) =1. In the limit 
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Fig. 5. Asymptotic value (for large Q 2 ) of the momentum fraction of the quarks 
inside the neutral pion (Eq. 26), as a function of the cut-off A. 

= 0, we obtain 



2 (x) (A = oo,to„. = 0) 



41n(2i//m 



4\n(2u/m 2 



(27) 



This means that in that regime, for sufficiently high Q 2 , the quarks behave 
as free particles, and the usual derivation based on the OPE holds [18]. 
The second case where this holds confirms this interpretation: if we impose 
m q < m n/2, the expressions we have given develop an infrared divergence, 
which corresponds to the case where both quarks emerging from the pion are 
on-shell and free. This divergence can be re-absorbed into the normalisation 
(25) of g, and the sum rule 2 (x) — 1 is again automatically satisfied at large 
Q 2 - 

However, in the physical pion case, it makes more sense to consider that one 
of the quarks remains off-shell: the imposition of a cut-off changes the sum 
rule value, as the fields can never be considered as free. Hence, because of the 
Goldstone nature of the pion, one expects that the momentum sum rule will 
take a smaller value than in the case of other hadrons. This may explain why 
fits that assume the same momentum fraction for valence quarks in protons 
and pions [14] do not seem to leave any room for sea quarks [19]. The results 
for 2 (x) are shown in Fig. 4. Again, the curves show a plateau at sufficiently 
large Q 2 , with a value depending on m q /A and m n /A. It is easy to show that 
this value is always smaller than 1. It is shown in Fig. 5 as a function of A. 

To fix the remaining parameters, we must then use our knowledge of con- 
stituent quarks, and of the pion. We choose conservative values m q = 300 MeV 
or 360 MeV, and a pion radius of the order of 0.25 fm, which corresponds to 



10 



m =0.300 GeV A=0.811GeV 
m '=0.300 GeV A=0.750 GeV 
in =0.360 GeV A=0.750GeV 



\ w 



\\\ 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Fig. 6. Structure function F2 for the neutral pion. On the left, Q 2 
values of the parameters m q and A are as indicated. On the right, 
A=0.75 GeV and the values of Q 2 are as indicated. 



2 GeV 2 and the 
0.3 GeV and 



A ~ 800 MeV. This choice of parameters gives us a momentum fraction 2 (x) 
between 0.55 and 0.65 (see Fig. 4), and corresponds to a coupling g with the 
plateau value of 3.8. It is remarkable that this value is very close to the one 
that guarantees in the NJL model, with the same constituent mass and the 
same cut-off, a unit value for the residue of the qq propagator at the pion pole 
and the correct value of the electric form factor at Q 2 = [17]. In other words, 
this corresponds to the value needed for the pion to appear in the NJL model 
as constituted of a quark and an antiquarkQ 

Let us finally examine the properties of the distribution v(x), or equiva- 
lently, of the function F2. Some of our results are summarized in Fig. 6, for 
Q 2 = 2 GeV 2 . The most striking feature is the vanishing of this function for 
x larger than some value x max . This is again due to the fact that, because the 
quarks are not free in this model, the actual value of their mass does mat- 
ter. It can clearly be seen that this effect originates from kinematical cuts: 
in the chiral limit and for large v, the condition s > 4m 2 is equivalent to 



1 



2m- 



in the absence of cut-off. 



In the case of a finite cut-off (A 2 -C Q 2 ), conditions (13-14) lower the value of 



71 



to 



X ■/■ 



1 

A 2 ' 



(2* 



This means that, for a finite Q 2 , when x is large, which corresponds to a small 
value of s, there is no way to put the cut quarks on their mass shell: this 
requires at least an energy of 4 m 2 . For given Q 2 and A, the available energy 

2 One should note that, although the numerical values are compatible, the cut-off 
has different physical origins in the two approaches: it corresponds to the maxi- 
mum internal momentum of the pion in our approach, while it is a parameter for 
regularizing ultraviolet divergences in the NJL model. 
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is increasing with decreasing x. As a result, putting the cut quarks on-shell 
will be easier for small than for large x. Therefore, the x-distribution (Fi) is 
expected to be enhanced on the low x side, leading to a momentum fraction 
smaller than unity. The vanishing at large x is not obtained in similar works, 
in particular in the one of Ref. [7]. In this reference, the Bjorken limit is taken 
first and the kinematical constraint (Eq. 13) is not applied . This procedure 
is certainly not correct for evaluating cross-sections at finite Q 2 . Except for 
this vanishing, our results basically agree with those of Ref. [7] . Therefore we 
will not present results for the DGLAP evolution of our structure function F 2 , 
which we expect to be compatible with the existing data. 



4 Discussion and conclusion 

We have discussed the simplest model allowing to relate virtual photon-pion 
forward elastic scattering to quark distributions. In this model, the imaginary 
part of the forward elastic scattering does not show any divergence. However, 
this quantity cannot provide the quark distributions readily since the numer- 
ical importance of the so-called crossed diagrams precludes the existence of 
such a relationship. The introduction of a cut-off for the relative momentum 
of the quarks inside the pion allows such an interpretation: crossed diagrams 
then appear as higher twists. The introduction of the cut-off does not allow 
to fulfill the momentum sum rule (2(x) = 1) at infinite Q 2 because, in the 
case of the pion, constituent quarks can never be considered as free. We have 
mentioned above how the kinematical constraint, allowing quarks to be put 
on their mass shell, leads to the reduction of the momentum fraction. This 
can easily be seen from our results for the chiral limit (Eq. 27). 

We motivated the cut-off as a manifestation of the pion size. The value 
A -1 = (0.75 GeV)^ 1 is close to the hard core rms radius of the chiral bag 
model, 0.35 fm [20]. 

The cut-off has been imposed on the relative quark momentum. This procedure 
is at variance with the double-subtraction Pauli-Villars procedure proposed in 
Ref. [9]. We have considered such a procedure in our model, but it produces 
discontinuities in the structure function as well as a negative value for this 
quantity in some regions of x. This procedure is also different from with the 
other ones introduced in similar works based on NJL models, where they are 
part of the necessary regularization of these models to avoid divergences. In 
addition, there is no need in our approach to consider additional diagrams with 
local pion-pion-quark-quark interactions. Yet, the pion-quark-quark coupling 
constant turns out to be the same in our approach and in NJL models. This 
may not be too surprising, as in both cases this coupling is determined by the 
requirement that the pion appears as made of two constituent quarks. 
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Our main conclusion is that pions are different from other hadrons, in the 
sense that the quark momentum fraction should be smaller, and that higher 
twist terms disappear for Q 2 ~ 2 GeV 2 . 



This work has been performed in the frame of the ESOP collaboration (Euro- 
pean Union contract N° HPRN-CT-2000-00130). We thank Dr. M. Diehl and 
Dr. P. Guichon for their useful comments. 
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